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ABSTRACT 


The  first  phase  of  a  terrain-clutter  study  has  been  completed  at 
NRL,  utilizing  the  four- frequency  radar  system  installed  in  a  WV2- 
128324  (Super  Constellation)  aircraft.  This  system  is  capable  of  trans¬ 
mitting  and  receiving  four  frequencies  (P,  L,  C,  and  X)  consecutively 
with  a  choice  of  horizontal,  vertical,  and  alternating  polarizations.  The 
returns  are  gated  at  a  fixed  range  and  analyzed  in  a  digital  computer. 
Absolute  clutter  measurements  have  been  taken  over  both  Canadian  and 
American  sites  at  angles  of  incidence  of  10  degrees  or  less.  Two  pa¬ 
rameters  have  been  measured.  The  first  is  the  normalized  radar  cross 
section  in  terms  of  the  10,  50  (median),  and  90  percentile  values  of  its 
distribution,  while  the  second  is  the  fluctuation  power  spectrum  of  the 
coherent  radar  return  signal.  Two  types  of  statistical  processes  have 
been  determined  from  the  cross-section  data;  one  follows  a  Rayleigh 
distribution  between  10  and  90  percentile  values,  and  one  departs  radi¬ 
cally  from  this  distribution.  The  energy  in  the  fluctuation  power  spec¬ 
trum  has  been  estimated  in  frequency  bands  whose  spacing  is  determined 
by  the  transmitted  frequency,  e.  g.  ,  2  cps  at  P  band  through  10  cps  at 
X  band.  The  root  mean  square  deviations  of  the  estimates  indicate  like¬ 
wise  a  Rayleigh  behavior  for  certain  bands  and  a  radical  departure  from 
this  behavior.  In  this  case,  the  radical  departure  is  identified  as  a 
Ricean  (signal  plus  noise)  distribution  produced  by  strongly  reflecting 
point  targets  in  the  illuminated  area.  A  statistical  model  for  terrain 
clutter  is  proposed  which  will  enable  the  phenomenon  disclosed  by  the 
study  to  be  simulated  based  upon  the  optical  analysis  of  Synthetic 
Aperture  Radar  Maps, 
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NRL  TERRAIN  CLUTTER  STUDY 
PHASE  I 


INTRODUCTION 

The  Naval  Research  Laboratory  is  currently  conducting  a  clutter  measurement  pro¬ 
gram  for  the  Research  and  Technology  Division,  Systems  Engineering  Group,  Wright- 
Patterson  AFB,  Ohio.  This  report  is  a  summary  of  results  and  data  accumulated  during 
the  first  phase  of  the  program,  covering  18  months  of  effort,  and  is  offered  in  lieu  of 
bimonthly  reports  four,  five,  and  six.  During  the  course  of  the  program  to  date,  three 
major  reports  (1, 2,  3)  were  produced  which  contained  large  quantities  of  data  and  anal¬ 
yses  pertinent  to  them,  which  were  circulated  to  workers  in  the  field. 

The  purpose  of  the  present  work  as  contained  in  the  MIPR  AS- 5-64  is  to  "contribute 
technical  data  to  establish  clutter  models  . . .  that  can  be  used  for  theoretical  evaluation 
of  various  rejection  techniques  ..."  To  achieve  this  purpose,  data  were  collected  by 
personnel  in  the  NRL  instrumented  aircraft  WV- 128324  (Fig.  1).  The  initial  objective  of 
the  study  was  to  acquire  data  on  four  frequencies  and  two  polarizations,  at  low  angles  of 
Incidence  (less  than  10  degrees)  over  both  American  and  Canadian  terrain.  The  data 
collected  were  to  be  processed  to  determine  both  the  amplitude  distribution  and  power 
spectra  of  the  clutter  return.  Due  to  operational  difficulty  in  acquiring  spectral  data, 
early  flights  yielded  only  partial  success;  however,  by  July  1965  most  of  these  difficul¬ 
ties  were  overcome,  and  the  complete  set  of  signals  was  recorded.  Table  1  gives  a 
listing  of  the  flights  on  which  the  majority  of  the  data  in  this  report  were  taken,  with 
comments  on  the  nature  of  the  data  acquired. 


Fig.  1  -  Radar-instrum.’nted  Supor-Constellr  ion  NRL  Aircraft 
used  in  terrain-clutter  studies 


The  major  portion  of  the  instrumentation  carried  by  the  NRL  aircraft  is  the  four- 
frequency  (4F)  radar  system,  a  research  radar  which  is  capable  of  transmitting  coher¬ 
ently  an  eight-pulse  sequence  of  the  form  P„,  X„,  L„,  C„.  P,.  X,,  L, ,  C,,  where 
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Table  1 

Flight  Summary 


Date 

Flight  Path 

Elevation 

Angle 

(degrees) 

Polarization 

Frequency 

Objort 

Sept.  8,  64 

Goose  Bay,  Labrador  to 

Winnipeg 

4-12-9 

VV,  HH,  VH.  HV 

X  band 

Cross  Section 

Sept  12.  64 

Nebraska,  Wyoming.  Utah 

5-7 

VV.  HH,  VH,  HV 

X  band 

Cross  Section 

Sept.  13,  84 

Milford,  Utah  to  SO. 

Arizona 

3-5 

VV.  HH.  VH.  HV 

X  band 

Cross  Section 

Sept.  18.  64 

Dec.  10.  64 

Dyersburg,  Tenn.  to 

West  Va. 

Eastern  Shore,  Md 

3— 1/2— 7 

4 

VV,  HH,  VH,  HV 

VV.  HH 

X  band 

X  and  C 
band 

Cross  Section 

Spectra 

■July  2,  65 

Iowa.  Ill  .  Ind.  ,  Ohio 
(Denver  flight! 

8 

VV,  VH.  VH.  HV 

P.  L.  C, 

X  bands 

Cross  Section 
and 

Spectra 

the  subscripts  refer  to  the  transmitted  polarization  The  repetition  rate  of  the  sequence 
is  variable,  with  a  maximum  rate  of  788  per  second,  so  that,  for  example,  the  prf  as¬ 
sociated  with  the  P„  transmission  is  788  pps. 

A  simplified  blor  k  diagram  of  the  receiving  system  is  shown  in  Fig.  2.  For  each 
pulse  transmission,  the  horizontal  and  vertical  components  of  the  backscatter  are  separ¬ 
ately  sensed  by  the  antenna  system,  converted  coherently  to  i-f,  and  fed  to  lin-log  and 
linear  i-f  amplifier  pairs  which  are  designed  to  preserve  amplitude  and  phase  of  the  en¬ 
velope  of  the  return  respectively.  In  the  logarithmic  strip  the  dynamic  range  is  com¬ 
pressed  for  easy  digitalization,  incoherently  detected  (square  law),  and  used  as  output  to 
the  range -gating  and  digitalization  circuitry.  Simultaneously,  in  the  linear  strip,  the 
signal  is  amplified  with  minimum  phase  distortion,  ha-d  limited  to  eliminate  amplitude 
variations,  coherently  detected,  and  is  used  as  output  o  the  range-gating  and  digitali¬ 
zation  circuitry.  The  other  polarization  component  is  treated  in  similar  fashion.  Some 
question  has  arisen  concerning  the  effect  of  the  hard-limiting  process  on  the  clutter 


Fig.  2  -  Simplified  block  diagram  of  four -frequency  radar 
receiving  system 
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spectrum;  this  will  be  treated  in  a  later  section  of  the  report.  For  the  next  pulse  trans¬ 
mission,  the  frequency-conversion  unit  and  the  stable  signal  source  are  switched  to 
provide  the  correct  signals  and  circuitry  to  convert  the  new  frequency/polarization  trans¬ 
mission  to  i-f,  and  the  same  lin-log  strip  system  is  used  for  amplification  and  detection. 
The  common  i-f  system  not  only  eliminates  receiver  errors  but  also  simplifies  the 
calibration  procedure.  A  detailed  account  of  the  4F  radar  system  will  be  found  in 
Appendix  A. 

A  nomenclature  has  been  evolved  to  describe  the  various  outputs  of  the  system.  The 
letters  P,  L,  C,  and  X  designate  the  transmissions  on  frequencies  of  428,  1228,  4455, 
and  8910  Me  respectively.  The  polarized  returns  will  be  symbolized  by  letter  combina¬ 
tions  of  the  form  HH,  VH,  etc.  ,  where  the  first  letter  denotes  the  transmitted  polariza¬ 
tion  and  the  second  the  received.  The  frequency -band  letters  will  be  used  for  the  ampli¬ 
tude  of  the  return  with  polarization  subscripts,  while  the  phase  of  the  returns  will  be 
designated  by  with  polarization  subscripts.  It  should  be  noted  that  for  the  eight -pulse 
sequence  previously  indicated,  there  are  32  recorded  outputs  per  range  gate,  producing 
a  data  rate  of  25  kilowords  per  second. 

The  transmitting  and  receiving  antenna  system  consists  of  four  dual-polarization 
antennas.  These  antennas  are  mounted  in  pairs  back  to  back.  The  X-  and  C-band  anten¬ 
nas  are  mounted  side  by  side,  and  the  beams  are  colinear  and  have  the  same  beam- 
width,  so  that  they' illuminate  the  same  area.  L-  and  P -band  antennas  are  mounted  180 
degrees  in  azimuth  from  X-  and  C-band  antennas,  and  are  likewise  colinear  but  have  dif¬ 
ferent  beamwidths,  so  that  the  area  illuminated  by  the  P-band  antenna  includes  the  area 
illuminated  by  the  L-band  antenna.  The  antenna  pairs,  while  locked  in  azimuth,  are 
separately  trained  in  elevation  from  0  degree  (horizontal)  to  90  degrees  (vertical).  The 
entire  array  may  be  scanned  through  315  degrees  of  azimuth,  or  any  portion  thereof. 

Two  views  of  this  system  are  shown  in  Fig.  3. 

Two  modes  of  data  collection  were  employed  in  the  flights  listed  in  Table  1.  The 
first  is  a  searchlight  mode,  in  which  the  antennas  are  held  at  fixed  azimuth  and  elevation 
angles.  These  are  held  relatively  constant  by  stabilization  referenced  to  vertical  and 
magnetic  north.  The  second  is  a  scanning  mode,  in  which  the  antenna  system  is  sector 
scanned  through  a  prescribed  angle  continually  for  the  duration  of  the  run.  In  both  cases, 
the  aircraft  is  attempting  to  maintain  a  constant  heading,  altitude,  and  velocity  over  the 
desired  flight  path.  The  searchlight  runs  are  used  primarily  to  acquire  spectral  data 
and  amplitude  variations  along  a  strip.  The  scanning  runs  allow  collection  of  the  ampli¬ 
tude  variation  statistics  over  a  larger  area  to  determine  wavelength  and  angle -of  - 
incidence  dependence.  In  order  to  describe  the  nature  of  the  terrain  along  the  flight 
path,  maps  have  been  provided  which  indicate  the  geographical  area  in  which  the  data 
were  taken  and  the  flight  path  as  determined  from  the  aircraft  instruments  and  visual 
observations.  These  maps  are  shown  in  Figs.  4  through  8.  The  flight  path  is  shown 
as  a  vector  labeled  with  the  run  number  (octal  notation).  The  vector  lies  along  the  ground 
track.  The  date  and  run  number  specify  a  unique  flight  path,  while  time  along  the  path 
is  specified  by  the  relative  time  clock  in  the  recording  system. 

The  remainder  of  this  report  will  deal  primarily  with  the  data  acquired  and  will  at¬ 
tempt  to  provide  sufficient  information  to  enable  the  reader  to  utilize  the  data,  within  its 
limitations.  In  addition,  certain  generalizations  will  be  made  concerning  the  nature  of 
clutter  and  the  clutter  process.  Since  the  primary  objective  of  this  program  is  to  imple¬ 
ment  a  clutter  model  for  radar  evaluation,  a  statistical  model  is  proposed  which,  it  is 
fek,  will  include  both  the  geographic  and  reflective  properties  of  the  terrain  in  a  manner 
suitable  for  computer  implementation. 


b.  X  and  C  dishe  s 


3  -  Four-frequency  antenna  system  installed 
in  lower  radome  of  aircraft 
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*»•  99*  98*  97*  96'  64*  63'  82*  6l*  SO*  79*  76*  77*  72*  7l*  70*  69*  68* 


Fig.  4  -  Flight  path,  Sept.  8,  1964,  Canada 


TERRAIN  CLUTTER 
Clutter  Analysis 

Clutter  Is  generally  considered,  from  the  standpoint  of  mathematical  tractability,to 
be  the  energy  scattered  back  to  a  radar  from  a  random  array  of  targets  of  roughly  equal 
cross  section.  The  scattei  ers  are  further  assumed  to  be  uniformly  di&tribated  over  the 
illuminated  area  and  along  the  flight  path.  The  signal  returned  is  the  result  of  the  vector 
addition  (random  walk)  of  the  contributions  of  each  scatterer  illuminated  by  the  radar. 

The  number  of  scatterers  is  assumed  sufficiently  large  to  ensure  convergence  of  the  sta¬ 
tistics.  This  set  of  assumptions  will  be  termed  the  "sandpaper  clutter  model”  for 
brevity.  Fluctuations  in  the  amplitude  of  the  clutter  signal  arise  a=  the  result  of  several 
mechanisms:  (a)  the  relative  motion  of  the  scattereis  with  respect  to  each  other  (inter¬ 
nal  motion);  (b)  variations  in  the  phase  of  the  returns  from  the  various  scatterers  with 
resp>.  t  to  the  radar  because  of  antenna  scanning  motion;  (c)  similar  variations  pro¬ 
duced  by  platform  motion  (aircraft  velocity);  (d)  phase  shifts  induced  by  error  motions, 
i.  e. ,  pitch,  roll,  and  yaw  of  airframe,  velocity  variations,  etc.  The  sandpaper  model 
has  been  analyzed  by  several  scientists  and  the  effects  of  each  of  the  above  mechanisms 
assessed.  Barlow  (4)  has  curve  fit  a  Gaussian  spectral  shape  to  internal- motion  phase 
variations  by  observing  the  fluctuations  of  a  searchlighting  fixed  radar.  The  variance  of 
the  distribution  is  given  for  various  situations,  wind-blown  leaves,  waves,  etc. 

Emerson  (5)  among  others  has  analyzed  the  fluctuation  spectra  produced  by  antenna  and 
aircraft  nut  ions  and  has  developed  Gaussian-shaped  spectra  for  their  phase  variations. 
The  antenna  is  assumed  to  be  Gaussian  in  shape-,  with  a  standard  deviation  equal  to  the 
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lh  the  above,  rc(0  is  the  unit  clutter  power  spectrum  whose  variance  is  <rc;  ain,  <ja  and 
apm  are  the  standard  deviations  of  the  internal  motion,  scanning,  and  platform  motion 
spectra  respectively;  v  is  the  standard  deviation  of  particle  velocities  in  the  clutter  cell; 
</,' ,  the  scan  rate;  u  the  aircraft  velocity;  a,  the  antenna  aperture;  the  antenna  angle 
measured  from  ground  track;  and  A,  the  wavelength.  Therefore,  the  normalized  clutter 
correlation  function  of  lag  r  becomes 

<Pr(r)  -  e  "  (770.  T>2  .  (3) 

Having  the  spectral  distribution  and  the  correlation  function  of  the  clutter  process,  it  is 
only  necessary  to  determine  the  total  distribution  of  the  cluster  amplitudes  to  describe 
the  process  completely.  The  sandpaper  model  is  a  form  of  the  two-dimensional  random- 
walk  problem,  which  yields  a  Rayleigh  distribution  for  the  ctutter  envelope  amplitude,  r, 


v 

PO  )  -  _  rxp 

rr  2 


(4) 


where  \  is  the  standard  deviation  of  the  amplitude.  From  the  radar  equation,  the 
clutter  cross  section  s  is  proportional  to  r2,  so  that 


(5) 
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Fig.  7  -  Flight  path,  Sept.  16,  1964,  Tennessee  to  West  Virginia 


Fig  8  -  Flight  path,  July  2,  1965,  Iowa  to  Ohio 
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where  the  nonrandom  radar  parameters  have  been  eliminated  by  normalization.  The 
radar-system  parameters  and  flight  geometry,  combined  with  the  above  description  of 
the  clutter  process,  enable  the  analyst,  at  least  in  theory,  to  predict  the  performance  of 
conventional  radars.  The  performance  characteristics  obtained  are  strongly  influenced 
by  the  values  used  for  the  variance  of  the  internal  motion  spectrum  and  the  value  used 
for  the  terrain  cross  section,  s.  Unfortunately,  these  parameters  vary  widely  over  the 
surface  of  the  earth,  and  exact  measurements  are  difficult  to  obtain.  For  radar  sys¬ 
tems  in  which  the  platform  motion  spectrum  has  a  small  variance,  the  internal  motion 
spectrum  has  a  strong  influence  on  performance;  however,  for  most  airborne  systems 
the  platform  motion  spectrum  is  sufficiently  large  to  dominate  the  clutter  spectrum,  so 
that  the  internal-motion  spectrum  has  little  influence  on  performance.  In  both  cases, 
the  spatial  variation  of  the  terrain  cross  section  is  of  primary  importance. 

The  specification  of  terrain  cross  section  is  complicated,  since  its  value  depends 
upon  the  roughness  of  the  terrain,  angle  of  incidence,  wavelength  of  the  incident  radia¬ 
tion  and  its  polarization,  and  the  size  of  the  illummated  patch.  The  cross  section  is  also 
a  function  of  the  azimuth  at  which  the  data  are  taken.  This  variation  is  handled  in  these 
data  hy  eiiher  a  tacit  assumption  of  azimuthal  symmetry,  where  a  fixed  angle  is  used, 
or  by  operationally  scanning  the  antenna  of  the  radar  during  data  collection  to  integrate 
azimuthal  effects.  The  following  sections  will  deal  with  measured  values  of  the  cross- 
section  distribution  and  the  clutter  power  spectrum,  both  to  enable  implementation  of  the 
Sandpaper  clutter  model  and  to  provide  values  for  analyses  where  the  assumptions  of  the 
model,  i.  e. ,  spectral  shape,  cross-section  distribution,  stationarity,  etc.,  are  suffi¬ 
ciently  imprecise  to  cause  difficulty. 


Clutter  Cross  Section 

The  basic  parameter  measured  in  the  amplitude  channel  of  the  four -frequency  radar 
is  the  logarithm  of  the  square  of  the  instantaneous  amplitude  of  the  clutter  return,  quan¬ 
tized  to  seven-bit  precision  (128  levels).  By  suitable  calibration,  described  in  Appen¬ 
dix  B,  an  absolute  value  of  the  terrain  cross  section  can  be  achieved.  In  processing  the 
data,  the  logarithmic  form  is  preserved  and  logarithmic  axes  are  used.  With  these  axes 
the  amplitude  and  the  cross  section  become  identical  and  may  be  used  interchangeably, 
a  point  which  has  caused  some  confusion  in  the  past.  The  value  of  the  cross  section  is 
given  by  the  radar  equation  as 


a 


(4’r  )3fl4rr 


\2G2(0,  VO  PT 


(6) 


where  the  radar  parameters  have  the  usual  significance.  The  constant  k  is  known  from 
the  system  parameters,  and  geometry,  and  ( Pr/PT )  is  determined  by  the  calibration  pro¬ 
cedure.  Tho  value  of  a,  thus  obtained,  is  dependent  upon  the  area  illuminated  by  the 
radar.  To  eliminate  this  dependence,  it  is  standard  to  compute  the  normalized  cross 
section,  defined  as 


•" 0 


(V 


where 


! 


co.s 


for  small  depression  angles. 
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■  in  v' 


for  large  depression  angles. 


(8) 


The  raw  dzta  which  are  recorded  In  flight  are  processed  in  a  small,  general-purpose 
digital  computer  (Raytheon  250)  in  the  following  manner.  A  relative  scale  is  established 
by  utilizing  the  fact  that  due  to  quantizat.on  of  the  signal  into  128  levels,  it  is  possible  to 
determine  the  cumulative  distribution  function  by  noting  the  number  of  signals  in  a  given 
sample  time  which  exceed  each  level.  The  scale  is  calibrated  by  adding  the  appropriate 
values  (in  decibels)  of  I ,  the  illuminated  area,  sample  size,  and  k  to  each  level,  thus 
yielding  the  absolute  value  of  The  level  which  is  exceeded  by  50  percent  of  the 
pulses  in  a  given  sample  time  is  the  median  value  of  the  distribution,  This  param¬ 
eter  was  chosen  as  characteristic  of  the  distribution  rather  than  the  mean  to  which  it 
is  related  In  a  Rayleigh  distribution  by 


.  1.6  JS. 


(9) 


since  it  is  less  sensitive  to  errors  which  are  concentrated  at  the  high  or  low  ends  cf  the 
distribution.  The  errors  at  the  high-signal  end  of  the  distribution  are  attributable  to  re¬ 
ceiver  saturation,  while  those  at  the  low -signal  end  are  attributable  to  averaging  in  the 
finite-width  range  gate,  and  receiver  noise.  To  minimize  the  gating  error,  a  30-nanosec 
gate  width  is  normally  used. 

For  each  time  sample  in  a  given  data  run,  the  entire  distribution  is  computed.  Two 
forms  of  presentation  of  the  cross-section  data  will  be  found  in  this  report.  For  runs  in 
which  the  antenna  is  scanned  at  constant  elevation,  the  entire  distribution  of  cross  sec¬ 
tions  is  given;  that  is,  the  total  ensemble  of  range -gate  samples  over  the  entire  run  are 
distributed  and  calibrated.  The  data  are  presented  in  tabular  form  with  selected  per¬ 
centile  values.  Normally,  the  100  percentile  is  the  lowest  received  pulse,  while  the  0 
percentile  is  the  value  of  the  largest  pulse  in  the  distribution.  Occasionally,  the  lower 
end  of  the  distribution  will  be  in  the  receiver  noise.  The  tables  will  th?n  commence  with 
the  percentile  level  at  the  minimum  detectable  signal. 

A  sampling  of  the  data  taken  in  this  program  may  be  found  in  Tables  2  through  6. 
Table  2  lists  the  operating  conditions  for  the  September  1964  flights,  while  Table  3  lists 
the  cumulative  distributions  for  the  various  runs.  A  similar  format  is  used  for  the 
Denver  data.  H  is  valuable  to  determine  from  these  data  the  validity  of  the  assumed 
Rayleigh  envelope  distribution.  A  convenient  representation  of  the  closeness  of  fit  of  the 
clutter  distribution  to  the  assumed  distribution  is  the  plot  of  data  points  on  Rayleigh  pa¬ 
per.  Figures  9,  10,  and  11  are  examples  of  this  representation.  The  ordinate  is  the 
value  in  decibels  of  the  amplitude  level,  while  the  abscissa  is  the  percentage  of  the  sam¬ 
ples  which  exceed  that  level  in  a  fixed  time.  The  sample  size  is  1000  pulses,  which  at 
788  pps  is  1.  27  sec  of  flight  time.  Figure  9a  shows  the  nature  of  the  fit  which  may  be  ob¬ 
tained.  A  Rayleigh  distribution  plotted  in  this  fashion  should  have  a  slope  equal  to  that  of 
the  theoretical  line  drawn  on  the  figure.  The  data  shown  are  roughly  parallel  to  the  theo¬ 
retical  line  through  the  middle  portion  of  the  plot  and  depart  at  the  high  and  low  ends  for 
reasons  previously  stated.  Figure  9c  shews  a  radical  departure  from  the  assumed  dis¬ 
tribution.  A  major  attempt  will  be  made  in  the  course  of  this  study  to  determine  the 
effect  on  a  generalized  radar  of  this  departure  from  Rayleigh  upon  the  probability  of  de¬ 
tection  and  false -alarm  rate.  The  remaining  figures  show  similar  plots  for  four  frequen¬ 
cies  and  two  polarizations. 
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Date 

Sept.  8,  64 


Sept.  12,  64 


Sept.  13,  64 


Sept.  16,  64 


Table  2 

Operating  Conditions,  Sept.  1964  Flights;  Antenna 
Scanning,  Transmitted  Polarization  V,  H 


Run 

Pulse  Width 
(  sec) 

PRF 

(PPS) 

Elev. 

Angle 

(degrees) 

Azimuth 

Angle 

(degrees) 

Sample 

Size 

55 

0.  25 

394 

5-1/2 

30  L  -  90  R 

42,800 

56 

4-1/2 

130  L  -  OR 

48,700 

57 

4-1/2 

60  L  -  60  R 

33,  400 

60 

5 

180  L  -  45  R 

19,500 

61 

5 

60  L  -  60  R 

32,500 

67 

9 

150  L  -  10  R 

24,700 

70 

4-1/2 

130  L  -  20  R 

27,700 

71 

8 

90  L  -  60  R 

46,  500 

103 

5-1/2 

70  L  -  80  R 

36,700 

104 

4-1/2 

70  L  -  80  R 

40,  800 

110 

0.  25 

788 

5 

70  L  -  70  R 

57, 500 

111 

7 

70  L  -  70  R 

107,400 

112 

7 

70  L  -  70  R 

46,  900 

113 

6 

70  L  -  70  R 

81,500 

114 

5 

70  L  -  70  R 

66,  200 

116 

5 

70  L  -  70  R 

82,700 

117 

5 

70  L  -  70  R 

83,  100 

120 

5 

70  L  -  70  R 

41,000 

161 

0.5 

788 

3 

70  L  -  70  R 

61,200 

162 

3 

70  L  -  70  R 

71,700 

163 

3 

70  L  -  70  R 

37,700 

165 

2 

70  L  -  70  R 

46,400 

166 

2-1/2 

70  L  -  70  R 

42,700 

167 

5 

70  L  -  70  R 

32,  400 

170 

5 

70  L  -  70  R 

83, 300 

64 

1.  0 

394 

5-1/2 

±50 

16, 600 

65 

7 

±50 

11,700 

66 

6-1/2 

110  L  -  50  R 

35,500 

67 

5 

110  L  -  50  R 

26,  500 

70 

5 

110  L  -  50  R 

26,800 

71 

4 

110  L  -  50  R 

17,800 

72 

3-1/2 

100  L  -  40  R 

30,700 

73 

3-1/2 

90  L  -  40  R 

21,500 

74 

3-1/2 

90  L  -  40  R 

41,500 

12 
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Table  3 

Cumulative  Distribution  of  Normalized  Cross  Section,  X  Band 
Sept.  8,  1964 
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Table  3  (cont'd) 
Sept.  12,  1964 


(Table  continues) 
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Table  3  (cont'd) 
Sept.  13,  1964 


(i  coutinues) 
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Table  3  (cont'd) 
Sept.  16.  1964 


Percentiles 


Run 

zatlon 

00 

75 

50 

25 

10 

0 

67 

VV 

-32.  50 

-25.  50 

-21.00 

-15.  25 

-10.00 

HH 

-24.  50 

-22.  50 

-20.  00 

-  5.50 

VH 

-24.  50 

-20.00 

HV 

-28.  00 

-20.  00 

70 

W 

-25.  50 

-21.00 

-18.00 

-11.50 

HH 

-25. 00 

-22.  50 

-20  50 

-15  00 

VH 

-25.  00 

-20.  50 

HV 

-29.  00 

-21.  00 

71 

VV 

-27.  50 

-22.  75 

-14.50 

HH 

-23.  50 

-21.  50 

-17.00 

VH 

HV 

-26.  50 

73 

W 

-26.  25 

-22.  50 

-15.  00 

HH 

-24.  50 

-22.  25 

-20.  50 

-16.  00 

VH 

1 

-23.  50 

-20.  00 

HV 

-22.  50 

74 

VV 

-29.  00 

-23. 00 

-15.  00 

HH 

-23.00 

-21.00 

-16.  00 

VH 

-24.  00 

-20.  50 

HV 

-21.75 

Table  4 

Operating  Conditions,  July  2,  1965  Flight;  Pulse  Width 
0.  5  «sec,  PRF  788  pps,  Elevation  Angle  8  Degrees 


Run 

Azim. 

Angle 

(degrees) 

Antenna 

Scanning 

Trans. 

Pol. 

Sample 

Size 

1 

60  L-60  R 

Yes 

V,  H 

138, 000 

2 

2  R 

No 

V.H 

93,000 

3 

25  R 

No 

V 

187,000 

4 

20  R 

No 

H 

40,  000 

6 

60  L-60  R 

Yes 

V,  H 

104,000 

7 

5  L 

No 

V,  H 

94, 000 

25  R 

No 

V 

115,000 

«W5*V-. 
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Table  5 

Cumulative  Distribution  of  Normalized  Cross  Section 
July  2,  1965 

Run  1 


Frequency  and 
Polarization 

Percentiles 

90 

75 

50 

25 

10 

0 

Xvv 

-33.  5 

-30 

-25.5 

-22.5 

-20 

-14 

Cvv 

-37 

-32 

-27.5 

-24.5 

-22 

-15.5 

Lvv 

-44.5 

-39 

•33.  5 

-30 

-26.  5 

-18.  5 

Pvv 

-50.5 

-44.5 

-39 

-35 

-32 

-22.5 

-34 

-30 

-25.5 

-22 

-19.  5 

-14.5 

C(|M 

-39 

-35 

-30 

-26 

-23 

-15.5 

Ltm 

-41.5 

-38.5 

-33 

-30 

-26 

-14 

P 

till 

-43.5 

-38.5 

-34 

-30 

_ 

-26.  5 

-14 

Run  2 


Frequency  and 
Polarization 

Percentiles 

90 

75 

50 

25 

...  ..  ....  ..  ,_| 

10 

0 

X,v 

-34.5 

-31 

-26 

-23.5 

-21 

-14.5 

Cvv 

-36.5 

-31.5 

-26.  5 

-24 

-21.5 

-15.5 

Lvv 

-42 

-36.5 

-32.  5 

-27 

-23 

-13 

Pvv 

-50.5 

-44 

-39 

-35 

-30.  5 

-21.5 

*H„ 

-33.5 

-29.  5 

-25.5 

-22.  5 

-19.  5 

-12.5 

C  |l  H 

-39 

-36 

-30,  5 

-27 

-24 

-17 

L||h 

-41 

-37 

-31.5 

-27 

-22 

P 

*  MM 

— 

-37 

3, 

L~ 

-25 

-22 

-16.  5 

i  TabU-  c ontinuc s) 
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Table  5  (rout'd) 

Run  3 


Run  4 


Frequency  and 
Polarization 


CVv 

L\  v 

Pw 


90 


75 


Percentiles 
50  25 


10 


V  Polarization  Not  Transmitted 


x„„ 

-35.5 

-31 

-27 

-24 

-21.  5 

-14 

c„„ 

-39 

-37 

-31.  5 

-28 

-24.  5 

-15.5 

“Mil 

*42 

-40 

-35.  5 

-31.  5 

-28.  5 

-15.  5 

P 

-44  j 

-39 

-32 

-27.  5 

-24.  5 

-  7.5 

X„v  j 

-42 

-38 

i 

-34 

-31.  * 

-29 

-22.5 

c 

H  V 

-44 

-41.5  j 

-37 

-33.  5 

-30.  5 

-23 

L„v 

-48.  5 

-46 

-42.  5 

-38.  5 

-35.  5 

-26.5 

P,n- 

- 

! 

-51 

-26.  5 

-43 

-32 

17 


oncy  and 
rization 

Percentiles 

90 

75 

50 

25 

10 

0 

x,, 

-34 

-30.5 

-25.  5 

-23 

-20 

-14 

C,v 

-37.5 

-32 

-27.  5 

-24.  5 

-21.5 

-15.5 

Lv  v 

-43 

-37.  5 

-33 

-29 

-26.  5 

-18.5 

P,v 

- 

-50.5 

-43 

-38 

-34 

-24 

X„„ 

c„„ 

H  Polarization  Not  Transmitted 

J-MII 

P„H 

X  vh 

-39.5 

-36.  5 

-33 

-29 

-26.  5 

-21 

C  VH 

- 

-39 

-37.  5 

-33 

-30 

-24 

"\|l 

-42.  5 

-41.  5 

-39.  5 

-37 

-29.  5 

P.H 

- 

-51 

-47 

-41 

-32.5 

(Tabl.  continues) 
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Table  5  (cont'd) 
Run  6 


Frequency  and 
Polarization 

Percentiles 

90 

75 

50 

25 

10 

0 

xvv 

-34.5 

-31 

-26 

-23 

-20 

-11.5 

Cvv 

-37 

-31.5 

-27 

-24 

-20.5 

-10 

Lvv 

-41 

-36.5 

-32.5 

-27.5 

-24.5 

-15 

Pvv 

-44 

-39.5 

-34.5 

-30 

-27 

-21 

X„H 

-33.5 

-29 

-24.5 

-21 

-16.  5 

-  3.5 

^HK 

-38.5 

-34.5 

-27.5 

-23 

-18.5 

-  3.5 

Lllll 

-42 

-38.5 

-34 

-29.5 

-25.5 

-12.5 

P 

mi 

1 

-47 

-41 

-36 

-32 

-27.5 

-18.5 

Run  7 


Frequency  and 
Polarization 

Percentiles 

90 

75 

50 

25 

10 

0 

Xvv 

-36 

-32 

-28.  5 

-25 

-22 

-14.5 

Cvv 

-37 

-32.5 

-28 

-24.5 

-22 

-15.5 

Lvv 

-40.5 

-36 

-32 

-27.5 

-26 

-18.5 

Pvv 

-45.5 

-41.  5 

-37.  5 

-34 

-31.5 

-23 

-35 

-31 

-27 

-24 

-21.5 

-15.5 

C  hh 

-39 

-36.5 

-30 

-27 

-24 

-17 

k>IH 

-41.5 

-37 

-33 

-28.  5 

-26 

-18.5 

p 

*»IH 

-42 

-38 

-34 

-31 

-28 

-22.5 

(Table  continues) 
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Table  5  (cont'd) 


Run  10 


Frequency  and 
Polarization 

|  Percentiles 

90 

75 

50 

25 

10 

0 

xvv 

-38 

-34 

-29 

-25 

-22.  5 

-16.5 

c 

v  vv 

-41 

-35.5 

-30.  5 

-25.5 

-23.5 

-17 

Lvv 

-40.5 

-36.5 

-32.5 

-27.5 

-26 

-18.5 

Pvv 

-44.5 

-39 

-34.  5 

-30.5 

-26.  5 

-21 

X,I„ 

c 

^  HH 

H  Polarization  Not  Transmitted 

P 

*  HH 

Xv„ 

*40 

-38.5 

-35 

-30.5 

-27.  5 

-21 

Cv„ 

- 

-39 

-"8 

-33.  5 

-30 

-23  5 

Lv„ 

-42.5 

-42 

-40.  5 

-36.  5 

-34 

-25 

PVH 

- 

-53 

-47.5 

-42 

-39 

-32.  5 

Table  6 

Operating  Conditions,  Dec.  10,  1964  Flight 


Run 

Pulsewidth 

(<sec) 

PRF 

(PPS) 

Elev. 

Angle 

(degrees) 

Azim. 

Angle 

(degrees) 

Antenna 

Scanning 

Trans. 

Pol. 

Sample 

Size 

0.  25 

1576 

4 

20 

V,  H 

33,096 

Using  the  July  1965  data,  the  wavelength  dependence  of  the  terrain  clutter  was  de¬ 
termined.  Two  presentations  of  the  dependence  are  demonstrated;  the  first,  Fig.  12a, 
has  the  range  of  values  of  the  normalized  ciutter  cross  section  as  a  function  of  the  wave¬ 
length.  Figure  12b  gives  the  ratio  of  the  clutter  cross  sections  as  a  function  of  the 
wavelength  ratio.  Six  such  ratios  may  be  obtained,  and  the  range  of  values  is  shown. 

The  solid  lino  in  each  figure  is  the  locus  of  values  which  would  be  obtained  for  an  in¬ 
verse  first  power  of  the  wavelength  dependence.  The  data  used  for  these  graphs  are  the 
median  values  of  clutter  cro3S  section  obtained  from  the  runs  taken  on  this  day.  The 
polarization  was  chosen  to  be  vertical  transmit  and  receive,  although  horizontal 
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Normalised  cross  section  (DB) 
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PERCENTAGE  Of  TIME  T*  UU.UES  ARE  lESS  THAN  THE  ORDINATE  IIOO  PI 
9999  99  90  60  40  20  10  5  2  I  05  01  005  0  01 


a.  Run  122,  elevation  angle  5  degrees 

9999  99  90  60  40  20  10  5  2  I  05  Oi  005  001 


PERCENTAGE  of  -ME  T*  VALUES  EXCEED  the  ORDINATE  (PI 
c.  Run  125,  elevation  angle  6  degrees 

Fig.  9  -  Measured  clutter  distributions 
compared  with  assumed  distributions 
(Rayleigh  slope);  Sept.  12,  1964,  VV 
polarization,  X  band 
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PE«a*ruf  of  tint  Ttc  me  less  thm  t*  omsart  noo  n 


a.  Polarization  VV,  X  band 


e.  Polarization  VV,  L  band 


PE«atT«£  of  r  ut  tk  values  eaceeo  the  qpC'Vats  ,pt 
g.  Polarization  VV,  P  band 


lEBCE N««  Of  T*  T*  VALUES  ME  LESS  tKAd  THE  OEWKTE  000  PI 


b.  Polarization  HH,  X  band 


i.  Polarization  HH,  L  band 


P£#tfNTA&£  Of  Tint  Inf  VALUES  £*C£E0  Tm£  OPClNATf  (P| 

h.  Polarization  HH,  P  band 


Fig.  11  -  Measured  clutter  distributions  compared  with  assumed 
distributions  {Rayleigh  slope);  July  2,  1965,  run  7 
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a.  Normalized  clutter  cross  section  as  a  function  of  wavelength 


x/p 


C/P 

X/L 


C/L 

L/P 

X/C 


b.  Ratio  of  clutter  cross  sections  as 
a  function  of  wavelength  ratio 


Fig.  12  -  Wavelength  dependence  of  terrain  clutter 

transmit  and  receive  show  little  difference  in  dependence  from  vertical.  This  result  is 
consistent  both  with  previous  NRL  measurements  and  also  those  taken  in  this  wavelength 
range  by  Ohio  State  University.  Much  of  the  previous  data  are  presented  by  L  Katz  of 
Applied  Physics  Laboratory  of  Johns  Hopkins  University  in  a  paper  entitled  ’Wavelength 
Dependence  of  the  Radar  Reflectivity  of  Earth  and  Moon,”  Vol  71,  No.  2,  Journal  of 
Geophysical  Research,  Jan  15,  1966.  Katz  has  shown  that  the  radar  cross  section  on 
the  average  is  a  decreasing  function  of  wavelength  to  approximately  the  first  power, 
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with  exceptions  whose  dependence  varies  from  +2nd  to  -6th  power  of  the  wavelength.  The 
positive  values  of  the  exponent  are  found  in  the  dependence  of  snow -covered  surfaces  at 
all  depression  angles  and  of  the  sea  surface  above  73  degrees  of  depression. 

R  is  well  known  that  the  clutter  cross  section  is  strongly  influenced  by  the  terrain. 
Two  types  of  variations  may  be  defined:  the  relatively  rapid  fluctuations  produced  by  the 
change  in  phase  of  the  scattered  energy  produced  by  internal  or  aircraft  motion,  and  the 
relatively  slow  variation  produced  by  gradual  changing  of  the  terrain  character.  These 
slow  variations  constrain  the  sample  time  over  which  a  distribution  may  be  considered 
stationary.  An  example  of  both  such  fluctuations  is  shown  in  Fig.  13.  This  run  was 
taken  over  Pisgah  Crater,  California.  In  order  to  specify  the  nature  of  the  scattering 
terrain,  photographic  imagery  was  used  both  down  looking  and  along  the  bore  site  of  the 
antenna.  This  picture  was  combined  with  USGS  aerial  survey  maps  of  the  area  to  pro¬ 
duce  a  flight  path  superimposed  on  the  terrain  synchronized  with  the  time  track  of  the 
recorded  data.  The  data  show  graphically  the  need  for  including  both  the  fast  fluctua¬ 
tions  and  slow  fluctuations  of  the  clutter  signals  in  any  model  which  is  to  be  representa¬ 
tive  of  the  terrain-clutter  return. 


Clutter  Spectra 

The  apparent  velocity  of  the  clutter  patch  illuminated  by  the  radar  is  given  by: 

*'  va  cos  d0  cos  c.;,  T  A'..-)  (10) 

where  va  is  the  aircraft  velocity;  and  are  the  elevation  and  azimuth  angles  of  the 
beam  center,  and  A0  is  the  angular  position  of  scatterers  in  the  patch  with  respect  to 
beam  center.  For  pulsed  radars,  the  variation  in  elevation  angle  over  the  patch  may  be 
neglected.  The  motion  of  the  patch  produces  a  Doppler  shift,/,  in  the  carrier  given  by: 

A/  -  2ty  \  cos  (<£„  +  At/’)  (11) 

where  «•,  is  the  ground  speed  of  the  aircraft,  and  \  the  transmitted  wavelength.  The 
variation  in  the  Doppler  shift  across  the  patch  gives  rise  to  a  spectrum  centered  at: 


whose  width  is  a  function  of  A ■/-■ .  For  a  nominal  ground  speed  of  300  ft/sec,  the  spectrum 
centered  at  X  band  has  a  maximum  value  of  5454  cps,  while  at  C,  L,  and  P  bands,  the 
center  is  located  at  2727  cps,  909  cps,  and  248  cps  respectively.  The  maximum  prf's 
available  in  the  4F  radar  system  are  1576  pps  (single  polarization)  or  1463  pps  (alter¬ 
nating  polarization).  If  the  radar  pulse  is  considered  as  sampling  the  Doppler  return, 
then  ambiguity  results,  since  the  sampling  rate  is  less  than  the  spectrum  center  fre¬ 
quency.  Doppler  or  velocity  ambiguity  is  a  possibility  in  three  of  the  four  frequencies. 

The  P-band  spectrum  is  unambiguous.  As  a  result  of  the  ambiguity  in  X.  C,  and  L 
(at  forward  angles),  spectra  will  exist  at  centers,, •„  ,  given  by: 

-■  U  *  “I'r  (13) 

where  /,  is  the  pulse -repetition  frequency.  Assuming  a  value  of  5454  cps  for  the  X-band 
Doppler  center  and  (r  equal  to  1576  pps,  representative  X-band  Doppler  spectra  will  be 
centered  at  726  cps,  850  cps,  2302  cps.  and  at  infinitely  many  other  locations.  It  is  char¬ 
acteristic  of  the  sampling  process  that  the  maximum  energy  is  located  in  the  spectrum 
centered  between  zero  and  one-half  prf.  R  is  also  characteristic  that  only  one  spectrum 
will  exist  in  this  range.  The  spectrum  centered  at  726  cps  differs  from  the  sampled  spec¬ 
trum  in  two  respects:  a  constant  amplitude  factor  and  a  constant  frequency  shift.  Since 
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TIME  (SECONDS) 

Fig.  13  -  Pisgah  Crater  terrain  cross  section,  Sept  13,  I  964,  run  14. 
elevation  angle  30  degrees,  polarization  VV,  X  band 


absolute  values  of  the  spectral  amplitudes  must  be  determined  by  calibration,  the  consta.it- 
amplitude  factor  will  be  absorbed  in  the  process,  while  the  original  frequency  axis  can  be 
reconstructed  from  a  knowledge  of  the  aircraft  velocity  and  geometry,  so  that  neither  of 
the  above-mentioned  characteristics  is  a  severe  liability.  The  spectrum  centered  at  726 
cps  is  termed  an  "alias"  of  the  original  spectrum,  since  only  the  frequency  label  of  the 
spectrum  has  changed  and  not  the  relative  shape.  The  spectra  presented  in  this  report, 
since  they  involve  C-,  X-,  and  L-band  transmitted  frequencies,  are  the  aliases  of  the 
original  spectra.  The  P-band  spectra  are  unique. 

If  the  width  of  the  spectra  is  such  that  portions  extend  beyond  the  limits,  1/2  prf  or 
prf.  a  further  imaging  of  the  spectra  will  result,  and  an  alias  of  the  energy  above  one- 
half  prf  will  appear  at  an  equal  frequency  below  prf.  This  phenomenon  is  termed  "fold- 
over"  and  destroys  the  usefulness  of  the  spectrum.  To  avoid  foldover,  the  alias  must  be 
positioned  roughly  midway  between  zero  and  one -half  prf,  since  foldover  also  occurs  at 
zero  frequency.  The  positioning  of  the  spectrum  is  achieved  in  the  field  by  the  adjust¬ 
ment  of  v0,  the  antenna  azimuth,  which  varies  the  apparent  velocity  of  the  clutter  patch 
arid  thus  the  frequency  of  its  Doppler  alias.  It  has  been  possible  with  X  and  C  bands  to 
find  an  optimum  azimuth  so  that  foldover  was  avoided  on  both  frequencies.  The  addition 
of  L  band  poses  a  more  difficult  task,  since  foldover  must  be  avoided  on  three  separate 
Doppler  aliases  simultaneously. 

Since  the  spectrum  position  is  a  function  of  the  position  angles  of  the  antenna  beam 
and  the  aircraft  velocity,  movements  of  the  antenna  or  changes  in  velocity  will  displace 
the  spectra  along  the  frequency  axis.  Averages  of  the  spectra  are  necessary,  so  that 
considerable  editing  of  the  data  is  required  to  find  portions  of  the  flight  path  where  the 
variations  of  the  aircraft  speed  and  antenna  motion  are  a  minimum.  The  criterion  pres¬ 
ently  used  is  that  a  long  time  average  of  the  spectra  must  converge  in  shape  to  the 
antenna  pattern  as  shown  in  Appendix  C.  This  criterion  eliminates  much  data:  however, 
that  which  remains  is  susceptible  to  statistical  analysis. 
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Two  forms  of  presentation  of  spectra  are  shown  in  this  report.  The  data  taken  over 
the  Eastern  Shore,  Maryland  are  displayed  with  a  linear  ordinate,  while  that  shown  for 
Denver,  Colorado  are  displayed  with  a  logarithmic  ordinate.  The  purposes  of  these  pre¬ 
sentations  differ.  The  Eastern  Shore  data  shown  in  Figs.  14  through  17  are  a  series  of 
spectra  which  were  taken  to  evaluate  the  effect  of  terrain  on  the  antenna  spectrum.  *  If 
the  antenna  is  subdivided  into  smaller  beams  of  width  A<//o,  then  the  illuminated  area, 

Eq.  (8),  is  subdivided  into  patches  of  width  R  a  <paf  where  R  is  the  range  to  the  patch  and 
length  c  t/2.  Equation  (11)  defines  a  doppler  spread  of  width  A  /„  over  the  smaller  area, 
where  the  total  of  all  such  frequency  bands  is  the  total  main  beam  spectrum.  The  power 
Incident  on  the  clutter  patch  is  given  by  PTr-  V'u,  ■/<„),  where  (>„  is  the  elevation  and  'P,, 
the  azimuth  of  the  center  of  the  patch.  If  the  patch  has  a  normalized  cross  section,  o-0„, 
then  the  power  in  the  clutter  spectrum  is 


p,  --  pT  £  a°“  °2  A/°  ■  (14) 

a  *  0 

Since  the  factors  r;  and  A fa  are  time  invariant,  the  power  variations  along  the  path  are 
related  to  the  values  of  so,  the  reflectivity  of  "ribbons"  of  terrain  of  width  R  A*0  paral¬ 
lel  to  the  flight  path.  In  order  to  determine  the  distribution  of  the  so  the  phase -versus - 
time  data  were  analyzed  by  digital  Fourier  analysis,  which  is  described  in  Appendix  C, 
and  the  standard  deviation  of  the  spectral  estimates  computed.  These  are  shown  in 
Figs.  18  through  21.  The  solid  curve  in  each  figure  is  the  value  of  r, J((',  <.'<),  the  two-way 
power  pattern  of  the  antenna  in  normalized  units.  The  bars  are  twice  the  standard  devi¬ 
ation  of  the  spectral  estimates  in  length. 

The  spectra  and  the  composite  spectra  document  the  existence  of  two  radically  dif¬ 
ferent  types  of  distribution  existing  in  the  terrain  ribbons.  The  first  is  a  random  distri¬ 
bution  whose  2<j  value  is  approximately  equal  to  the  mean,  a  characteristic  possessed  by 
the  Rayleigh  distribution.  The  second  is  a  distribution  which  from  inspection  of  the  spec¬ 
tra  is  produced  by  strong,  sustained  signals  present  in  the  clutter  and  may  be  termed  a 
"signal  plus  noise  distribution"  in  the  sense  of  Rice  (6). 

To  proceed  further,  the  noncoherent  envelope  correlation  function  is  defined  (7)  as 


7.  j(t)  * 


_  2 

M/’> cos  27,foTdf  + 


relf')  s>"  277  fo^f 


(15) 


where  ff  r  is  given  in  Eq.  (1),  f0  is  the  spectrum  center,  and  r  the  tag  time.  The  nor¬ 
malized,  but  not  dc  suppressed,  incoherent  video  correlation  function  is  given  by 


<t>c  (  T)  r 


1  +  7.  2 (  T) 
2 


(16) 


^Frequency  in  these  and  all  subsequent  figures  is  measured  in  cycles  per  second. 
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Fig.  19  -  Standard  deviation  of  spectrum  estimates,  run  57,  Dec.  10,  1964,  time  19013, 
Eastern  Shore,  Md.  ,  X  band,  VV  polarization,  10  degrees  of  freedom,  10  cps  bandwidth 
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while  the  pulse -to-pulse  phase  shift  is  given  by 


/'„(T)  *  MB*1 


(17) 


The  pulse  -to-pulse  phase  shift,  ,  ,  produced  by  a  point  target  in  the  clutter  patch,  is 
given  by  a  similar  expression,  where  the  spectrum  of  the  target  is  inserted  in  Eq.  (15). 
Urkowitz  (8)  has  shown  that  the  combined  video  correlation  function  of  a  fixed  point  target 
of  amplitude  P  in  a  clutter  patch  whose  average  power  is  i '  is 


T)  r 


(  «  *  I)1  +  2  X  cos  4  /*(’  ) 


(s  4  |>>  .  2«  4  I 


(18) 


where  t  P1 '2i 2. 

The  above  expression  indicates  that  the  pulse -to-pulse  noncoherent  correlation  func¬ 
tion  is  strongly  affected  both  by  the  amplitude  and  position  of  the  strong  point  scatterer. 
The  coherent  correlation  function,  which  is  proportional  to  the  square  root  of  the  inco¬ 
herent  function,  is  similarly  affected.  Since,  as  will  be  noted  in  the  next  section,  clutter- 
cancellation  techniques  use  the  {xilse -to-pulse  correlation  to  reject  clutter,  the  knowledge 
of  both  the  Rayleigh  -distributed  clutter  patches  and  the  signal-plus -noise -distributed 
patches  is  extremely  important  to  the  prediction  of  radar  performance  in  clutter. 

The  spectra  shown  in  Figs.  22  through  29  are  spectra  taken  on  the  Denver  trip,  July 
1965.  The  logarithmic  axes  used  to  present  these  spectra  preserve  the  maximum  dynamic 
range  available  in  the  quantized  output  of  the  phase  detector.  The  major  interest  in  these 
spectra  is  their  uniqueness,  in  that  in  conjunction  with  the  cumulative  distributions  given 
in  Table  5  for  this  run,  they  present  a  nearly  complete  statistical  description  of  the  clutter 
process  over  this  terrain  as  a  function  of  wavelength  and  polarization.  In  addition,  the 
correlation  function  which  may  be  derived  from  these  spectra  enables  the  cancellation 
ratio  to  be  determined  over  a  wide  variety  of  circumstances. 

In  order  to  assist  other  analysts  in  the  use  of  this  data,  a  brief  description  of  the 
clutter  power  spectra  presented  in  this  section  is  included  here.  This  description  in  con¬ 
junction  with  Appendixes  A  and  C  should  yield  sufficient  information  to  allow  general  use 
of  the  information.  Table  6  gives  the  operating  conditions  for  the  Eastern  Shore  flight. 

In  the  spectra  shown  in  Figs.  14  through  17  an  analysis  window  width  of  5  cps  was  chosen. 
This  value  was  established  by  using  a  sample  length  of  0. 1  sec.  This  value  of  window 
width  was  chosen  primarily  to  allow  rapid  processing  of  the  data.  Specification  of  the 
resolution  in  a  more  exact  manner  requires  a  knowledge  of  the  spatial  dependence  of  the 
terrain  features  which  at  present  is  absent,  since  imagery  of  the  terrain  is  not  available. 
Doth  X-  and  C-band  returns  were  analyzed  using  the  same  width.  Estimates  of  the  spectra 
were  made  at  every  10  cps,  In  X  band  from  100  to  500  cps,  and  in  C  band  from  10  to  250 
cps.  A  record  of  21  seconds  length  was  available  over  which  the  antenna  position  was  rel¬ 
atively  stable.  This  time  was  determined  by  the  test  mentioned  previously.  It  was  de¬ 
cided  to  eliminate  the  rapid  fluctuations  produced  by  the  changing  phase  of  the  scattererc, 
to  average  over  ten  consecutive  spectra  of  0. 1-sec  duration  to  form  the  resulting  curves. 
This  method  also  has  the  advantage  of  making  each  spectrum  the  average  over  *he  illumi¬ 
nated  area  of  the  beam,  since  the  aircraft  travels  a  distance  approximately  equal  to  t  -2 
(0.25-  sec  pulse)  in  one  second.  Due  to  a  defect  in  the  computer  input  program,  a  gap 
of  approximately  one  third  of  a  second  exists  between  each  graph  pliXted.  The  spectra 
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Fig.  25  -  Clutter  power  spectra,  July  2 
L  band,  VV  polarization 


1965,  run  2 
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Fig-  28  -  Clutter  power  spectra,  July  2,  19b5.  run  2, 

X  band.  HH  polarization 
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which  are  averaged  to  produce  the  plots  are  contiguous,  however.  Each  plot  is  labeled 
with  the  run  number,  polarization  transmit -receive  combination,  number  of  degrees  of 
freedom  in  the  average  (number  of  individual  graphs),  the  window  width  in  cycles,  the  ref¬ 
erence  clock  time  on  the  tape  recorder,  the  date  of  the  run,  the  frequency  transmitted, 
and  the  approximate  geographic  area.  The  frequency  axis  is  the  image  (alias)  axis,  not 
the  true  doppler  axis.  The  ordinate  of  the  spectra  is  relative  energy  (amplitude  squared). 
The  Denver  flight  spectra  (Figs.  21  through  28)  have  window  widths  of  10  cps  for  X-  and 
C-band  spectra,  5  cps  for  L-band  spectra,  and  2  cps  for  P-band  spectra.  These  spectra 
are  likewise  averages  containing  40  degrees  of  freedom  in  the  X-  and  C-band  spectra, 

20  in  L  band,  and  eight  in  P  band.  These  correspond  to  segment  lengths  of  0.  05,  0.  05, 

0.  1  and  0.25  second  respectively  for  the  individual  spectra.  In  all,  42  seconds  of  data 
were  analyzed,  resulting  in  21  spectra  for  each  of  the  frequency-polarization  pairs.  The 
ordinate  on  these  curves  is  decibels  above  the  system  noise  level,  which  may  be  assumed 
to  be  between  85  and  90  dbm. 

The  pointing  angles  involved  were  elevation  4  degrees,  azimuth  22  degrees  for 
Eastern  Shore,  while  for  Denver  they  were  8  and  0  degrees.  The  aircraft  velocity  is  a 
nominal  300  ft/sec. 


Cancellation 

A  conventional  delay-line  canceller  (9)  utilizes  a  delay  line  of  length  /  equal  to  the 
interpulse  interval  and  a  subtraction  device  to  form  the  difference  between  successive 
pulses  (single  delay)  or  between  successive  pulse  differences  (double  delay).  If  is 
the  amplitude  of  the  first  such  pulse  and  K  the  amplitude  of  the  second,  then  the  average 
residue  after  cancellation  of  a  number  of  such  pairs  is 


h  1  (R„  -  ;t  ) l'  t  *  /;  Kj  '  -  h^k  (19) 

where  t  is  the  expectation  operator.  The  cancellation  ratio  of  a  delay-line  device  is 
defined  as 

clutter  pt'Wr.  in  A  H  ‘ 

- - - - - -  (20) 

i  tuner  out  f-  '  (  N  H  )  '  ' 

If  /  i  .V;  >  f  t  k  ■  )  h  (  k  :  i  .  then  combining  Eqs.  (19)  and  (20) 

i  i 

>  - - whrrr  :.  (•) -  (21) 

.’<  1  -  »-n  h  H  ') 

and  combining  with  Eq.  (16) 


<  -  - -  (22) 

i  -  1  •  i 

Expanding  (Ref.  7)  1  in  terms  of  its  moments,  and  assuming  symmetry  of  the 

spectrum  so  that  only  even  central  moments  need  bo  considered. 


i 


(23) 
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And  therefore  for  main-beam  clutter, 


For  a  single  delay-line  canceller  r  r  r. 


where  cr:  is  the  variance  of  the  clutter  spectrum  in  radians  per  second,  and  f,  is  the 
pulse -repetition  rate.  To  determine  the  effect  of  approximating  /.'2(-0  by  its  second 
moment  along,  a  calculation  was  made  using  a  cosine -squared  approximation  to  the  4F 
radar  X -band  antenna  with  the  feature  that  the  side  lobes  are  adjustable.  The  pattern 
was  of  the  form: 


i(-TTT)1  fi  ‘ 

so  that 


c  — 


(24) 


(25) 


r;J(0)  =  <;a  cos2  (45/2  0)  -4°  <  0  £  4° 
r;2(0)  -  c,0  sin2  (450)  -3°  1  0  -  -4° 

4°  £  0  1  8°.  (26) 


Where  r;0  is  a  normalizing  constant.  Since  the  antenna  azimuth  pattern  is  related  to  the 
spectrum  of  a  pulsed  radar  by 


c2(0)  J4<--  *■(()  df. 


(27) 


Eq.  (15)  was  solved  with  the  aid  of  Eq.  (11)  for  \  =  0. 11  ft,  u  =  300  ft/sec,  and  T  = 
1/1576  second  for  as  a  function  of  the  side-lobe  level  The  result  was  converted 

to  cancellation  using  Eq.  (22).  The  result  of  this  calculation  is  shown  in  Fig.  30.  For 
side-lobe  levels  greater  than  10  db  below  the  main  beam,  the  increase  in  cancellation 
is  almost  linear  with  decreasing  side  lobes.  Consequently,  care  must  be  utilized  in 
the  calculation  of  the  cancellation  ratio  from  the  main-beam  clutter  spectral  width 
alone. 


CONCLUSIONS 

The  NRL  clutter  study  has,  in  the  first  phase  of  the  work,  accomplished  its  primary 
objective  of  measuring  clutter  cross  sections  at  small  angles  of  incidence  over  various 
American  and  Canadian  terrains.  Data  have  been  collected  on  four  frequencies  and  two 
polarizations.  In  addition,  spectra  have  been  measured  at  the  various  frequency- 
polarization  combinations.  In  the  process  of  the  study,  a  software  library  of  special 
programs  has  been  developed  for  data  processing,  and  more  effective  calibration  pro¬ 
cedures  have  been  developed.  Certain  theoretical  studies  have  been  performed  both  to 
determine  the  accuracy  of  the  data  and  to  provide  pertinent  parameters  for  specific  radar 
problems.  No  new  theoretical  concepts  have  been  discovered;  however,  certain  param¬ 
eters  whose  nature  has  been  hitherto  assumed  have  been  measured  and  presented.  The 
data  presented  in  this  report  are  not  considered  a  complete  coverage  of  the  clutter  prob¬ 
lem.  Obviously,  more  variety  of  terrain-type  studies  are  required;  for  example,  such 
categories  as  wooded  areas,  plains,  suburban  areas,  urban  areas,  mountains,  and 
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Fig.  30  -  E£fect  of  side-lobe 
level  on  cancellation 


deserts.  Transition  areas,  plains  to  mountains,  land  to  sea,  desert  to  mountains,  etc. , 
also  have  interest.  In  each  of  these  instances,  both  distributions  and  spectra  are 
required. 

In  the  course  of  the  study,  several  pertienent  questions  have  arisen  which  require 
answers  if  radar  performance  is  to  be  estimated.  What,  for  example,  is  the  effect,  if 
any,  of  departure  of  the  terrain- clutter  distribution  from  the  Rayleigh  distribution,  and 
what  is  the  effect  of  the  strong  point  scatterers  on  the  various  radar  processing  schemes? 
In  general,  the  treatment  of  the  entire  area  of  the  terrain  effects  on  radar  perfor¬ 
mance  which,  to  date,  has  consisted  of  a  single  parameter  to  describe  clutter  power  and 
an  assumed  spectral  form  given  by  the  sandpaper  model,  which  generally  ignores  ter¬ 
rain  differences,  is  obviously  inadequate  to  define  the  behavior  of  the  more  sophisticated 
systems  presently  under  consideration.  To  fill  this  need  for  specific  answers,  an 
accurate  and  flexible  clutter  model  is  required  which  will  contain  not  only  the  terrain 
characteristics  but  also  the  flexibility  to  include  the  parameters  of  the  radar  to  be 
analyzed. 


In  this  regard,  it  is  necessary  to  define  two  general  types  of  mathematical  models, 
the  rigid  model  and  the  statistical  model.  The  rigid  model  consists  of  deterministic 
expressions  obtained  by  physical  analysis  or  by  empirical  curve  fitting  of  data  which, 
when  combined  with  certain  input  parameters,  yield  a  prediction  of  performance  in 
terms  of  a  set  of  output  parameters  of  assigned  value.  A  statistical  model  (dynamic 
model)  consists  of  a  series  of  probabilistic  expressions  derived  from  physical  analysis 
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or  from  real  data  which,  when  combined  with  the  input  parameters,  yield  a  statistical 
measure  of  the  performance,  and  includes  the  rigid  model  as  a  special  case.  That  is, 
the  averages  derived  from  the  statistics  should  yield  the  rigid-model  results,  assuming 
sufficient  averaging  time  is  provided.  Since  the  terrain  seen  by  an  airborne  radar 
returns  energy  which  fluctuates,  both  because  of  changes  in  the  relative  phasing  of 
scattered  signals  and  because  of  variations  in  terrain  reflectivity,  bcth  of  which  are 
essentially  statistical  quantities,  a  model  which  is  able  to  represent  these  quantities  as 
data  would  seem  superior.  The  implementation  of  such  a  model,  while  difficult,  is  not 
impossible. 

In  the  course  of  this  study,  the  concept  of  a  statistical  model  based  upon  synthetic  - 
aperture  mapping  has  evolved.  A  synthetic -aperture  map,  or  for  that  matter,  PPI  maps 
and  scanning-aperture  maps,  are  plots  of  terrain  reflectivity  (radar  cross  section  per 
resolution  cell)  as  a  function  of  the  range  and  distance  along  the  flight  path  (azimuth  for 
PPI).  The  photographic  medium  offers  high  storage  capacity  combined  with  extremely 
slight  dependence  upon  the  viewing  radar.  The  parameters  of  primary  interest  are  the 
viewing  angle  and  resolution  cell.  So  within  these  limits  they  are  characteristic  of  the 
terrain  alone.  Such  a  map  can  be  converted  to  data  suitable  for  computer  processing  in 
a  manner  which  preserves  not  only  the  relative  values  of  reflectivity  but  also  the  spatial 
distribution  by  equipment  already  procured  by  NRL.  The  four-frequency  radar  system  is 
being  equipped  with  recorders  to  allow  the  synthetic -aperture  technique  to  be  applied  to 
all  signals  supplied  by  this  system. 

Briefly,  the  radar  maps  will  be  quantized  into  areas  equal  to  the  resolution  cell  of 
the  map  by  use  of  a  microdensitometer  whose  output  will  be  digitalized  to  provide  both 
a  measure  of  cross  section  and  coordinates.  The  resulting  information  will  be  recorded 
on  magnetic  tape  for  computer  entry.  In  the  computer  a  model  antenna  will  be  moved 
over  the  terrain  data,  and  each  cross  section  will  be  assigned  a  random  phase  and 
summed  within  the  illuminated  area.  As  the  radar  proceeds  over  the  terrain,  the  illu¬ 
minated  area  will  move  over  the  model  terrain,  and  thus  a  clutter  signal  as  a  function 
of  time  will  be  generated. 

The  major  advantage  of  this  technique  is  the  close  association  between  the  clutter 
data  and  the  geographical  location  of  the  scatterers.  Correlations  can  be  achieved  which 
were  beyond  the  elaborate  program  attempted  with  the  Pisgah  crater  data  described 
earlier  in  this  report.  The  major  disadvantage  is  the  difficulty  involved  in  calibrating 
a  radar  map.  A  more  detailed  analysis  of  this  technique  will  be  made  in  the  course  of 
the  second  phase  of  the  study.  E  should  be  noted  that  the  addition  of  the  mapping  equip¬ 
ment  will  in  noway  hinder  the  data  collection  facilities  presently  available. 

In  summary,  the  NRL  clutter  study  has  (a)  cataloged  data  over  a  number  of  terrains 
concerning  the  radar  cross  section  and  fluctuation  spectra,  (b)  presented  analyses  which 
demonstrate  the  existence  of  at  least  two  processes  in  terrain  clutter,  the  random 
Rayleigh-like  clutter  and  the  signal-plus -noise -type  clutter,  (c)  proposed  a  statistical 
model  based  on  synthetic  aperture  mapping  which  will  enable  the  evaluation  of  the  inter¬ 
action  of  radar  systems  with  specific  terrain  features. 
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Appendix  A 

THE  FOUR-FREQUENCY  RADAR  SYSTEM 


The  four -frequency  radar  system  is  basically  four  different  pulsed  coherent  radars 
transmitting  at  P  band  (428  Me),  L  band  (1225  Me),  C  band  (4455  Me),  and  X  band 
(8910  Me).  The  transmitted  signals  are  generated  at  a  low  level  by  a  stable  signal 
source  which  utilizes  a  1. 15  Me  precision  oven -controlled  Manson  basic  oscillator, 
varactois,  and  transistors  to  provide  multiplication  and  signal  amplification  to  the  de¬ 
sired  frequency  ranges.  In  addition  to  the  transmission  frequencies,  a  37 -Me  reference 
frequency  is  provided  for  use  in  coherent  detection.  The  outputs  of  the  stable  sigaal 
source  (30  mw)  are  fed  to  their  respective  transmitters  and  amplified  to  the  final  output 
levels.  Two  types  of  coherent  transmitter  chains  are  used,  namely,  tetrodes  in  external 
coaxial  cavities  for  P  and  L  bands,  and  one  or  two  TWTs  (rf  pulsed)  driving  klystron 
power  amplifiers  for  C  and  X  bands.  At  P  and  L  bands,  in  order  to  obtain  adequate  band¬ 
width,  the  driver  stages  are  operated  cw  and  the  output  stage  is  pulsed  with  switch  tubes 
to  provide  the  desired  rise  times. 

Each  of  the  four  radars  is  designed  to  operate  with  two  antennas,  one  polarized 
horizontally  and  one  vertically,  either  separately  or  pulsed  alternately.  The  four  fre¬ 
quencies  can  be  transmitted  singly  or  pulsed  in  rapid  succession  to  provide  a  total  of 
eight  different  frequency-polarization  combinations.  The  basic  timing  for  the  system  is 
provided  by  the  synchronizer  and  ranging  unit.  It  supplies  modulator  triggers,  video 
pulse  drives,  sweep  triggers,  timing  marks,  pulse -repetition  frequencies,  pulse  widths, 
and  gate  widths.  The  reference  used  for  this  unit  is  an  82-kc  clock  which  is  counted 
down  to  the  desired  rates.  In  each  radar,  an  antenna  switch  controlled  by  the  synchro¬ 
nizer  switches  the  transmitter  output  pulses  alternately  between  the  duplexer  systems  of 
the  two  antennas,  horizontal  and  vertical,  associated  with  a  given  frequency  transmission. 

The  four -frequency  antenna  system  consists  of  two  pairs  of  antennas  located  back  to 
back  on  a  common  mount.  One  pair  contains  the  lower-frequency  antennas,  P  and  L, 
which  are  crossed  dipoles  with  a  common  ground  plane.  The  second  pair  contains  the 
C-  and  X-band  antennas,  which  have  separate  parabolas  facing  180  degrees  away  from 
the  dipole  arrays.  The  antennas  can  be  rotated  through  315  degrees  of  azimuth  and 
100  degrees  of  elevation,  and  are  stabilized  in  roll  and  pitch  with  reference  to  magnetic 
north  and  vertical.  They  are  capable  of  being  sector  scanned  at  30  degrees  per  second. 
The  antenna  is  mounted  in  a  radome  on  the  underside  of  the  aircraft. 

On  reception,  the  signals  from  the  horizontal  and  vertical  antennas  are  passed 
through  their  respective  duplexers  to  the  two  microwave  mixers  and  preamplifier  units, 
where  they  are  converted  to  i-f  frequency  (37  Me)  and  amplified.  Each  of  the  eight  (two 
per  frequency)  i-f  preamplifiers  has  a  conventional  Wallman  cascade  input  circuit,  fol¬ 
lowed  by  a  buffer  tube  and  a  beam -deflection  tube  switch  to  time  share  the  common 
horizontal  and  vertical  i-f  receiving  channels.  The  horizontal  or  vertical  i-f  receiving 
channels  are  identical,  and  each  consists  of  parallel  i-f  amplifiers,  one  limiting  and  one 
logarithmic.  The  logarithmic  amplifier  has  a  bandwidth  of  12  Me,  a  lin-log  transfer 
characteristic  linear  within  0.5  db  over  45  db  of  dynamic  range,  and  performance  inde¬ 
pendent  of  overall  radar  duty-cycle  variations  of  600  to  1.  The  limiting  amplifier  is 
designed  to  produce  minimum  phase  error  over  a  large  range  of  input -amplitude  varia¬ 
tion.  The  amplifier  uses  a  Z-2800  beam -deflection  tube  which  has  excellent  limiting 
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characteristics,  constant  input  and  output  characteristics,  and  insignificant  input- output 
coupling.  Measurements  on  this  amplifier  have  shown  an  output  variation  of  a  small 
fraction  of  a  decibel  for  70  decibels  of  input  range.  The  phase  variation  was  within 
±2.  5  degrees  for  60  db  of  input  range. 

Two  forms  of  second  detection  are  used.  The  logarithmic  amplifier  has  a  detector 
with  a  square- law  characteristic,  while  the  limiting  amplifier  has  a  linear  synchronous 
detector  using  the  37-Mc  signal  from  the  stable  signal  source  as  reference.  An  addi¬ 
tional  phase  detector  is  included  which  allows  the  phase  difference  between  the  horizontal 
and  vertical  signals  to  be  measured.  Consequently,  five  video  outputs  are  available  for 
recording  for  each  transmitted  frequency,  namely,  the  amplitudes  of  the  dim.  i  and  cross 
polarized  return,  their  phases  with  respect  to  the  37-Mc  reference,  and  their  phase  with 
respect  to  each  other:  these  outputs  are  capable  of  being  displayed  on  A  scopes  for  the 
purpose  of  selecting  a  range  at  which  the  signal  is  to  be  recorded.  Two  range  gates  may 
be  used:  the  main  range  gate  which  may  be  set  with  an  accuracy  of  several  yards  out  to 
maximum  range  and  an  auxiliary  gate  which  may  be  positioned  with  respect  to  the  first. 
Digital  readout  to  the  nearest  yard  is  provided  for  the  main  range  gate. 

Finally,  a  16-head,  one- inch  magnetic  tape  system  is  used  for  data  recording.  The 
maximum  data  rate  is  25,  600  seven-bit  binary  words  per  second.  Each  of  these  words 
is  the  amplitude  existing  instantaneously  at  the  detector  outputs  at  the  time  of  the  range 
gate.  These  amplitudes  are  digitalized  in  groups  of  four,  since  four  amplitudes  are 
available  for  each  pulse  transmitted.  In  addition  to  these  high-rate  signals,  ten  other 
slowly  varying  parameters  may  be  recorded,  one  every  second.  Tables  Al,  A2,  and  A3 
give  a  listing  of  the  characteristics  of  the  system  which  has  been  outlined  and  give  a 
graphic  presentation  of  the  extreme  flexibility  embodied  in  the  four- frequency  radar. 


Table  Al 

Principal  Parameters  of  Four-Frequency  Radar  System  -  Transmitters 


— 

Parameter 

Frequencies 

428  Me 

1225  Me 

4455  Me 

8910  Me 

Long-term  frequency 
stability  -  no  greater 
than 

1 

in  106  per  day 

1  part 

in  10°  per  day 

1  part 

in  10(>  per  day 

1  part 

in  10(l  per  day 

Short-term  frequency 
stability  (0.  01  sec) 

Within  5  parts 
in  10"  -  10s 

Within  5  parts 
in  109  -  10“ 

Within  5  parts 
in  109  -  10“ 

Within  5  parts 
in  10';  •  10“ 

Peak  power 

35- 40  kw 

35-  40  kw 

25  kw 

40  kw 

Average  power 

140  watts 

140  watts 

100  watts 

160  watts 

Peak  power 
special  mode 

10  watts 

10  watts 

10  watts 

10  watts 

In ter  pulse  phase 

stabilit'’ 

no  pre.itcr  than 

4 

4  ' 

5 

6 

Output  Tube 

7651  tetrode  (RCA) 

! 

7051  tetrode  (RCA)  j  SAC  -  290 

i  Klystron  (Sperry) 

V24C  Klystron 
(Variant 

_  i 
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Table  A2 

Principal  Features  of  Four-Frequency  Radar  System  - 
_ Signal  and  Display  Parameters _ 


Parameter 

Value 

Repetition  rates 

, 

Sulimultiples  of  81,  959  kc. 

100.45,  200.90,  301.35,  394.07,  512.30, 

602.71.  683.07,  788.16,  1463,  2926  cps. 

Pulse  widths 

0  1.  0.25,  0.5,  1,  2/;  sec 

Noise  figure 

9-  11  including  all  rf  losses 

I-  F  amplifier 

Center  frequency 

Bandwidth 
,  Type 

Quantity 

i 

_ 

37  Me 

10,  4,  2,  1,  0.  5  Me 

lln-log,  25  and  50  db  dynamic  range 

Limiting 

Two  -  One  for  horizontally  polarized  signals, 
one  for  vertically  polarized  signals. 

Range  gate 

Width 

Increment  readout 

0.  1.  0.  2,  0.  5,  1,  3,  5/isec 

1  yard 

Rant!®  marks 

Width 

Spacing 

1.0,  2,0,  5.0,  20^sec 

I,  5,  20  nautical  miles 

Videi)  system 

Consoles 

Displays 

CRT  size 

Total  3:  2  operators  and  1  master  control 
Dual-beam  A  scope,  B  scope  and  PPI  lor  each 
console.  All  flat  face  tubes. 

5  in. 

Signal  outputs  (ungated 
and  gated) 

i _ 

Vertical  polarization  signal  ampliiudc 

Horizontal  polarization  signal  amplitude 

Vertical '  horizontal  polarization  phase 

Horizontal  polarization,  transmitted  signal  phase 
Vertical  polarization  transmitted  signal  phase 

1  Maximum  Video  Amplitude  j  2  volts,  peak-positive 

Table  A3 

Project  Suspender  Antenna  System  Performance 


Parameter 

P-Band 

L-Band 

C-BBnd 

X-  Band 

Polarization 

Horiz. 

Vert. 

Horiz. 

Horiz. 

Vert. 

Horiz. 

Vert. 

Bandwidth  (Me) 

±  5 

±  5 

t  5 

±10 

±10 

±10 

±10 

Azimuth  beamwidth 
(degrees) 

12.  3 

12.  1 

5.  5 

5 

5 

5 

4.  7 

Elevation  beamwidth 
(decrees) 

40 

41 

13.  8 

13 

5 

5 

5.  3 

5 

Azimuth  minor  lobe  (db) 

14.  5 

14.  5 

13.  4 

14 

23.  2 

23.0 

23.  6 

23.  6 

Elevation  minor  lobe  (db) 

30 

26 

16 

14 

24.  5 

24 

23.  5 

24.  2 

Cross  polarization  (db) 

25 

28 

25 

>20 

>20 

>20 

>20 

Transmission  iine  loss  (db) 

0.  6 

0.  3 

0.  6 

2.0 

2.  0 

0.  6 

Antenna  gain  (db) 

17.  4 

17.  4 

25.9 

31.  4 

31.  4 

31. 2 

31.  2 

Voltage  standing- wave  ratio 

<1.  44:1 

<1.  38:1 

<1 . 35:1 

^  1.28:1 

<  1.  50 

<  1.  57 

cl.  30:1 

1.  12:1 

High  power 

OK 

OK 

OK 

OK 

Not  Measured 

OK 

OK 

Appendix  B 

SYSTEM  CALIBRATION  AND  SOURCES  OF  ERROR 


Two  measurements  are  required  for  an  absolute  system  calibration:  (a)  a  measure¬ 
ment  of  the  receiver  linearity  and  gain,  and(b)  an  absolute  power  measurement  indepen¬ 
dent  of  the  system.  The  receiver  gain  and  linearity  measurements  were  made  using  an 
external  signal  generator  and  variable  precision  attenuators.  In  addition,  a  re-entrant 
test  signal  generated  within  the  system  is  monitored  during  the  data-collection  runs  to 
provide  a  check  of  receiver  stability.  Calibrations  were  thus  obtained  for  each  trans¬ 
mitted  frequency  polarization  signal. 

Absolute  power  measurements  were  obtained  by  tracking  an  eight -inch -diameter 
aluminum  sphere  dropped  from  the  plane.  The  sphere  is  weighted  to  fall  at  a  constant 
angle  relative  to  the  antenna.  The  magnitude  of  the  sphere  return  is  recorded  and  the 
range  logged.  Applying  the  radar  equation  to  the  sphere  return  denoted  by  the  sub¬ 
scripts)  yields 


s  - 


6-2 


R* 


TS 


For  the  terrain  measurement,  the  radar  equation  becomes 


sT  - 


(4n)> 


Tt 


The  ratio  of  Eq.  (B2)  to  Eq.  (Bl)  yields 


Ks 


(Bl) 


(B2) 


(B3) 


where  the  peak  transmitted  power  is  the  same  for  both  measurements.  Since  all  other 
quantities  are  measured,  the  terrain  cross  section  becomes  a  function  of  the  returned 
signal.  l‘Ki. 

There  are  two  basic  errors  in  this  calibration  system:  (a)  receiver  stability,  and 
(b)  tracking  error  in  the  sphere  measurement.  Previous  data  have  shown  that  the  sphere 
measurement  has  an  error  of  approximately  0  to  2  db,  due  chiefly  to  the  variation  of  the 
position  of  the  sphere  within  the  beam.  This  error  is  minimized  by  measuring  several 
spheres  for  each  calibration  and  using  the  strongest  return  at  a  given  range.  The  re¬ 
ceiver  stability  is  monitored,  as  stated  previously,  by  recording  a  known  reference  level 
at  various  time  intervals.  The  reference  has  an  uncertainty  of  approximately  1  db  for 
all  frequencies  and  polarizations. 

Consequently,  for  any  given  value  of  cross  section  a  maximum  error  of  ±3  db  is 
possible.  However,  the  polarization  ratio  would  contain  only  the  stability  error,  since 
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the  tracking  error  is  common.  Frequency  ratios  of  X  to  C  band  or  L  to  P  band  would 
likewise  have  only  stability  errors,  since  the  sphere  is  illuminated  by  these  frequencies 
simultaneously. 

The  above  analysis  does  not  take  into  account  all  possible  sources  of  error,  such 
as  variations  in  transmitter  power,  antenna  isolation,  receiver  crosstalk,  etc.  These 
errors  are  negligible  except  where  specifically  mentioned.  Since  the  received  signal 
is  measured  with  respect  to  a  sphere,  and  since  all  frequencies  and  polarizations  are 
calibrated  with  the  sphere  data,  comparisons  between  any  two  components  will  contain 
errors  principally  due  to  the  two  factors  discussed  previously. 


Appendix  C 
DATA  PROCESSING 


In  order  to  determine  the  Doppler  (clutter)  spectrum  from  the  phase -versus-timo 
samples,  two  approaches  are  possible  -  the  indirect  approach,  in  which  the  autocorrela¬ 
tion  function  of  the  data  is  first  determined  and  the  power  spectrum  derived  from  its 
cosine  transform,  and  the  direct  approach,  in  which  the  complex  amplitude  function  of 
the  data  is  determined  and  the  power  spectrum  derived  from  its  square.  The  latter  ap¬ 
proach  was  chosen,  since  it  offers  more  flexibility  in  the  choice  of  analysis  bandwidth 
and  frequency  resolution  and  allows  effective  monitoring  of  the  spectral  position. 

For  data  which  are  periodically  sampled  at  a  rate  f, ,  the  pulse-v  .  ;>-tition  rate,  which 
produces  a  time  spacing  between  samples  of  .',7.  the  Fourier  coefficients  «<  t )  and  mm 
are  given  by: 


1  > 

o(f)  —  .V  (n  A  T  )l  In  A  7)  cos  2’'f(n"T) 

K  .0 

(Cl) 

1  v 

6(f)  r  -  £  \(/»  A  7  )  ft(n  '  T)  sin  2  ’tin  A  T  ) 

ft  m  u 

where  V  is  the  number  of  points  in  a  sample,  '  7i  is  the  series  of  phase  points,  and 
*<"  *  T)  Is  the  sampling  function.  Several  choices  of  sampling  function  are  available; 
however,  in  view  of  the  quantities  of  data  to  be  processed,  a  sampling  function  which 
optimizes  program  speed  was  chosen,  namely: 


I  (n  T  )  ”  I  n  _  N 

(C2) 

0  otherwise. 

The  analysis  window  established  by  this  function  is  given  by  its  Fourier  transform: 


0 


(C  3) 


from  which  the  power  transfer  function  becomes: 


I  lid)!2 


sin ‘  t  T m 


(  "I  7  1  ’ 


rm 


\  •  r 


(C4) 


where  l„  is  the  length  of  the  sample.  The  equivalent  width  of  the  filter  between  zeroes 
is  then: 

’  f  '  7 f  (C5) 
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This  relationship  between  the  window  width  (equivalent  analog  filter  width)  and  the  sample 
length  allows  the  resolution  of  the  spectrum  to  be  controlled  simply. 

The  estimate  of  the  energy  in  the  spectrum  at  frequency  /  is  given  by: 

•i Hf)  -  «2(.0  i  (C6) 

and  the  estimate  of  the  total  power  in  the  spectrum  for  a  sample  length  Tm  is  given  by: 

fr 

P(f.  rn)  ~  2  j  A2d)df  •  (C7) 

0 

The  proportionality  exists  since  the  spectrum  is  uncalibrated,  while  the  factor  2  Is  in¬ 
cluded  since  both  positive  and  negative  frequencies  contribute  to  the  integral.  The  total 
power  in  the  present  system  is  a  constant,  since  the  limiter  provides  a  constant -input 
amplitude  to  the  phase  detector,  while  the  phase  between  the  reference  oscillator  and 
the  signal  is  sufficiently  variable  to  average  to  a  constant  multiplying  factor.  This 
system  characteristic  may  be  used  to  check  the  accuracy  of  the  computation.  The  total 
power  P<f.  10  rm)  was  computed  for  each  of  the  spectra  of  the  X-band  return  (HH 
polarization)  of  run  57,  Eastern  Shore,  Maryland  and  for  ar  equal  number  of  spectra  of 
the  C-band  return  (VV  polarization)  of  the  same  run.  The  percentage  deviation  from  the 
mean  was  5  percent  for  the  X-band  power  and  9  percent  for  the  C-band  power.  These 
results  were  considered  quite  satisfactory,  since  only  the  portion  of  the  spectra  con¬ 
taining  clutter  energy  was  used  in  the  computation.  This  portion  amounted  to  approxi¬ 
mately  50  percent  of  the  ranges  from  0  to  1/2  f,  for  X-band  and  25  percent  of  the  range 
for  C-band.  The  inclusion  of  the  system  noise  included  in  the  remainder  of  the  range 
would  lessen  the  percentage  deviation  considerably. 

A  second  test  of  the  accuracy  of  the  data  is  based  on  the  limiting  form  of  the  clutter 
spectrum.  The  limiting  form  may  be  determined  for  the  case  of  an  antenna  whose  two- 
way  power  pattern  is  c2(0)  which  is  modulated  by  any  terrain  cross  section  which  is 
independent  of  0  in  the  sampling  interval  and  which  has  a  distribution  given  by: 

P(S)  r  1  e-s/s°  s  >  o.  (C8) 


The  limiting  form  of  the  clutter  spectrum  becomes: 
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where  the  angle  axis  of  the  power  pattern  has  been  converted  to  a  frequency  axis  by  the 
transformation: 
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Equation  (CIO)  states  that  the  limiti.ig  form  of  the  clutter  spectrum  should  be  similar  in 
shape  to  the  antenna  power  pattern.  This  relation  was  verified  by  averaging  consecutive 
spectra  and  comparing  the  average  with  the  antenna  power  pattern  converted  to  frequency 
abscissa.  The  power  pattern  was  fitted  to  the  data  by  multiplying  it  by  an  arbitrary  con¬ 
stant  and  adding  a  correction  for  system  noise.  The  closeness  of  the  fit  is  shown  by 
Figs.  Cl  through  C4. 

Considerable  discussion  has  arisen  in  the  past  few  months  concerning  the  effect  on 
the  clutter  spectrum  of  the  hard  limiting  which  precedes  phase  detection.  In  general, 
hard  limiting  produces  harmonic  distortion  of  the  limited  signal.  In  the  present  system, 
hard  limiting  occurs  at  i-f  and  is  followed  by  a  bandpass  filter  centered  at  the  i-f  (37  Me). 
The  phase  distortion  due  to  limiting  is  approximately  5  degrees  of  phase  over  the  whole 
amplitude  range.  That  is,  the  phase  of  an  input  signal  is  measured  to  this  accuracy  in¬ 
dependently  of  its  amplitude.  The  pertinent  question  (aside  from  the  effect  of  the  5  degree 
errors)  then  is  how  much  information  is  common  to  the  phase  spectrum,  the  amplitude 
spectrum,  or  a  mixture  of  the  two.  McGinn  and  Pike  (7)  give  all  of  the  relationships 
that  they  could  find  for  narrow -band  noise-type  signals  and  point  out  (page  67)  that  the 
second -probability  densities  of  amplitude  and  phase  are  functions  of  the  same  single 
variable: 
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which  is  very  nearly  (page  70)  the  autocorrelation  of  the  signal  an  plitudes.  The  phase 
and  amplitudes  are  well  correlated,  but  the  author  knows  of  no  analysis  of  what  types  of 
signal  fluctuations  are  not  common  to  the  two  spectra.  This  correlation  has  been  empir¬ 
ically  demonstrated  by  side -looking  mappers  which  are  usually  greatly  limited,  yet  pro¬ 
duce  terrain  maps  which  are  in  effect  plots  of  terrain  cro^s  section  versus  planar 
coordinates.  As  a  result  of  these  arguments,  the  spectn.  .  resented  in  this  report  are 
an  accurate  measure  of  the  phase  spectrum  and  highly  correlated  with  the  amplitude 
spectrum. 
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Fig.  CZ  -  Average  spectrum,  run  37,  Dec.  10,  1964,  time  19013, 
Eastern  Shore,  Md.  ,  X  band,  VV  polarization,  10  degrees  of 
freedom,  10  cps  bandwidth 
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Fig.  C3  -  Average  spectrum,  run  57,  Dec.  10,  1964,  time  19013, 
Eastern  Shore,  Md.  ,  C  band,  HH  polarization,  10  degrees  of 
freedom,  5  cps  bandwidth 
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Fig.  C4  -  Average  spectrum,  run  57,  Dec.  10,  1964,  time  19013, 
Eastern  Shore,  Md.  ,  C  band,  VV  polarization,  10  degrees  of 
freedom,  5  cps  bandwidth 
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Department  of  the  Air  Force 
(Systems  Engineering  Group,  SEAFW) 

_ _ _ Wright- Patterson  Air  Force  Base 

The  first  phase  of  a  terrain -clutter  study  has  been  completed  at  NRL,  utilizing 
the  four -frequency  radar  system  installed  in  a  WV2-128324  (Super  Constellation)  air¬ 
craft.  This  system  is  capable  of  transmitting  and  receiving  four  frequencies  (P,  L, 
C,  and  X)  consecutively  with  a  choice  of  horizontal,  vertical,  and  alternating  polariza¬ 
tions.  The  returns  are  gated  at  a  fixed  range  and  analyzed  in  a  digital  computer.  Ab¬ 
solute  clutter  measurements  have  been  taken  over  both  Canadian  and  American  sites 
at  angles  of  incidence  of  10  degrees  or  less.  Two  parameters  have  been  measured. 
The  first  is  the  normalized  radar  cross  section  in  terms  of  the  10,  50  (median),  and 
90  percentile  values  of  its  distribution,  while  the  second  is  the  fluctuation  power  spec¬ 
trum  of  the  coherent  radar  return  signal.  Two  types  of  statistical  processes  have 
been  determined  from  the  cross-section  data;  one  follows  a  Rayleigh  distribution  be¬ 
tween  10  and  90  percentile  values,  and  one  departs  radically  from  this  distribution. 
The  energy  in  the  fluctuation  power  spectrum  has  been  estimated  in  frequency  bands 
whose  spacing  is  determined  by  the  transmitted  frequency,  e.  g. ,  2  cps  at  P  band 
through  10  cps  at  X  band.  The  root  mean  square  deviations  of  the  estimates  indicate 
likewise  a  Rayleigh  behavior  for  certain  bands  and  a  radical  departure  from  this  be¬ 
havior.  In  this  case,  the  radical  departure  is  identified  as  a  Ricean  (signal  plus 
noise)  distribution  produced  by  strongly  reflecting  point  targets  in  the  illuminated 
area.  A  statistical  model  for  terrain  clutter  is  proposed  which  will  enable  the 

(Abstract  continues) 
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phenomenon  disclosed  by  the  study  to  be  simulated  based  upon  the  optical  analysis  of 
Synthetic  Aperture  Radar  Maps. 


